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ABSTRACT 

The  effect  of  a finite  number  of  items  being  screened  is 
evaluated,  both  in  the  case  where  all  parameters  of  a bivariate 
normal  distribution  are  known  and  where  all  piurameters  are  un- 
known. Some  illustrative  tables  are  included. 
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1 . INTRODUCTION 

A survey  of  prediction  Intervals  and  their  applications  was 
given  by  Hahn  and  Nelson  (1973) . More  recent  work  has  been  done 
by  Fertig  and  Mann  (1975),  (1977a)  and  (1977b).  In  the  present 
paper  we  are  Interested  only  in  normally  distributed  variables. 
This  paper  differs  from  previous  papers  on  prediction  intervals 
for  the  nomud.  distribution  in  that  we  have  avail2d}le  observations 
on  a variable,  X,  which  we  will  use  to  screen  items  for  Inclusion 
in  a set  of  items  which  are  required  to  meet  a specification  on  a 
correlated  variable,  Y.  That  is,  we  will  consider  an  item  to  be 
good  if  Y ^ U,  and  we  want  to  be  able  to  assiire  that  the  number  of 
items  meeting  this  specification  is  at  least  k in  the  group  found 
acceptable  by  screening  on  X.  The  problem  is  to  set  a limit  on  X 
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which  will  accomplish  this  with  a given  probability.  We  will 
treat  two  cases,  one  where  all  parameters  are  known,  and  one  where 
all  parameters  are  unknown. 

2.  CASE  WHERE  ALL  PARAMETERS  ARE  KNOWN 

Owen,  et  al.  (1975)  proposed  a method  of  using  a variable  X 
correlated  with  a variable  Y to  screen  items  so  that  the  propor- 
tion of  y U is  raised  from  y before  screening  to  5 after  screen- 
ing. The  variables  (X,  Y)  are  assuned  to  have  a joint  bivariate 

normal  distribution  with  means  (u  , u ) , respectively;  standard 

X y 

deviations,  (a  , a ) respectively  and  correlation  p which  we  will 
X y 

assume  is  positive.  CVen,  et  al.  (1975)  give  tables  of  a quantity 
0 so  that  if  all  itasts  are  accepted  for  which  X < u 4-  K„  a then 

X D X 

the  goal  of  raising  the  proportion  of  Y's  less  than  V from  y to  i5 
is  accomplished.  The  quantity  is  the  normal  deviate  that 
corresponds  to  a proportion  6 in  the  lower  tail  of  a univariate 
normal  distribution. 

Owen,  et  al.  (1975)  point  out  that  the  proportion  <5  is  an 
expectation  achieved  for  an  entire  normal  population.  If  only  a 
finite  number  m of  items  are  screened  then  the  number  of  items,  V, 
meeting  the  specification  Y U is  a random  variable  following  the 
binomial  distribution  with  parameters  (m,  0) . Hence,  if  we  wanted 
the  P{V  ^ t)  C,  then  we  need  to  solve  the  equation 


I n«^(l  - «)“"^  - C 

J-f  ^ 

where 

6 - P{Y  £ olx^  - ‘'x  ’‘®x' 

for  k.  Note  that  we  cm  solve  the  first  equation  for  5 and  obtain 


where  F,  . is  100(1  - C)%  upper  tail  percentage  point  of  the  F- 

<;,a,o 

distribution  with  a degrees  of  freedom  for  the  numerator  and  b 


degrees  of  freedom  for  the  denominator.  Hence,  after  computing 
the  5 from  this  formula  we  proceed  to  compute  k as  in  Owen,  et  al. 
(1975)  from 

P{Y  < u|x  < u + ka  } - 5. 

The  result  is  shown  in  the  accompanying  t^d}les  under  infinite  de- 
grees of  freedom  for  a sample  of  size  100,  and  some  selected  values  of  ? . 

Note  that  the  adjustments  can  be  made  as  in  Owen  et  al.  (1975) 
for  both  lower  specification  limits,  L,  and  negative  correlation. 
Hence,  the  procedures  and  teUsles  given  here  apply  to  any  case  where 
there  is  a one-sided  specification  limit  and  either  a positive  or 
negative  correlation.  The  steps  given  above  also  apply  to  two-sided 
specification  limits  down  to  the  last  expression  for  5.  We  are 
then  100C%  sure  that  at  least  t of  m future  values  of  Y will  be 
less  than  U when  items  are  selected  based  on  X < u + ka  . 

— X X 

3.  CASE  WHERE  ALL  PARAMETERS  ARE  UNKNOWN 

In  Owen  and  Su  (1977)  the  following  procedure  was  given  to 
take  care  of  the  case  where  (u  , U , P,  a , a , y)  are  unknown. 

X y X y 

1.  Take  a preliminary  saiq>le  (Xj^,  y^^) , (Xj,  y^) , ...,  (x^, 
y^)  and  estimate  p,  a^,  o^)  in  the  usual  way,  i.e. 

__  n _ n 

X “ I X /n,  y - ^ y./n, 
i-1  i.l  ^ 


- I (Xj^  - x)  V(n  - 1) , Sy  - I (y^  - y)V(n  - 1) , 

i=l  ^ i-1 


I (x^  - x)  (y  - y) 

i-1 

(n  - 1)  s s 
X y 


2.  Compute  a 100ti%  lower  confidence  limit  y*  for  Y “ P{Y^U} 

by  finding  the  noncentrality  parameter  K * for  a noncentral  t- 

Y 

variate,  T^,  that  satisfies 


whare  k - (U  - y)/Sy.  Since  la  the  univariate  normal  deviate 
corresponding  to  y*  in  the  lower  tall  we  can  easily  obtain  Y*  and 
we  have 


p{y  > Y*)  ■ n- 


3.  Compute  a 100n%  lower  confidence  limit  0*  for  p.  This 
may  be  done  using  David's  <19S4)  tables  or  approximately  from 


P*  ■ tanh  < arctanh  r - 


»n 


Obtain  the  value  of  k from  Owen  and  Haas  (1978)  by  enter- 

. = ■ Y*. 


ing  those  teibles  with  D.F.  ■ n - 1,  D 


5.  For  all  additional  items  accept  those  for  which  X ^ x 


♦ a 


/n. 

/ n 


X 

6.  We  can  then  be  100 (2n  -1)%  sura  that  at  least  1005%  of 
the  Y values  in  the  screened  population  are  less  than  U. 

Again,  if  there  is  only  a finite  ntxnber,  m,  in  the  screened 
population  then  the  number  of  items,  V,  meeting  the  specification 
Y U follows  a binomial  distribution  with  parameters  (m,  5) . 
Again  we  want  P{\'  it}"  C wtiere  t is  some  minimum  nxanber  of 
items  meeting  specifications  which  we  want  to  see  among  the  m 
items  which  have  been  screened,  and  C is  the  probability  of  our 
seeing  this  result. 

The  procedure  is  the  same  as  in  Section  2 above  except  that 
this  time 


5 • P{Y  i U|x  j<  X )t  Sjj) 


where  5 la  computed  from  the  F-distribution,  as  before. 

This  time  the  quantity  on  the  right  must  be  obtained  from  the 
normal  conditioned  on  t-dlstribution  instead  of  the  normal  condi- 
tioned on  normal  distribution  as  in  Section  2.  Computational 
algorithms  are  given  by  even  and  Haas  (1978)  for  this. 


4 . TABLES 


In  the  accompanying  tables  we  give  values  of  k as  defined  in 
Sections  2 and  3 for  ra  ■ 100;  y ■ 0.8,  0.9;  p - 0.90,  0.99, 

; - 0.90  and  0.99,  i.e.,  just  a few  illustrative  values  since  the 
tables  would  bo  very  massive  to  cover  any  reasonable  range  of  uses 


TABLE  I 


Values  of  k When  m ■ 100  and  ; = 0.90 


n - 

1-30 

n - 1 

- infinity 

0 ■ 0.90 

0 - 0.99 

0 - 0.90 

p - 0.99 

X 

X 

.4 

40 

1.176 

1.176 

1.135 

1.136 

60 

.232 

.285 

-.227 

-.278 

80 

-.311 

-.073 

-.303 

-.071 

90 

— 

— 

-.620 

-.237 

.8 

80 

1.695 

1.722 

1.626 

1.642 

90 

.924 

1.110 

.909 

1.084 

.9 

90 

1.899 

1.968 

1.826 

1.868 

TABLE  II 


Values  of  k 

When  m ■ 100  and 

; - .99 

n 

-1-30 

n - 1 - 

infinity 

p - .90 

p - .99 

p - .90 

p - .99 

1 

X 

.40 

40 

.800 

.805 

.777 

.781 

60 

.092 

.176 

.090 

.172 

80 

-.427 

-.129 

-.377 

-.126 

90 

— 

— 

-.743 

-.257 

.80 

80 

1.341 

1.412 

1.303 

1.357 

90 

.739 

.999 

.729 

.987 

.90 

90 

1.520 

1.669 

1.486 

1.606 

5.  CONCLUSION 


Me  can  be  100(2n  - 2 + ?)  % sure  that  at  least  1 of  m future 
observations  on  the  variate  Y will  be  below  U if  the  m observations 
on  X all  have  been  screened  so  that 


X < X + k s / ~ . 

— X n 

The  quantity  k can  be  read  from  the  accompanying  table  for  some 
special  cases.  However,  in  most  instances  it  will  be  necessary  to 
compute  k using  the  methods  of  Owen  and  Haas  (1978)  for  the  case 
where  the  normal  conditioned  on  t probability  is  set  equal  to  5 as 
defined  above. 
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